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ABSTRACT 



The results of seven different tests of a single stage 
axial turbine with converging-diverging stator nozzle are 
reported. From measurements in a special test rig, losses 
occurring in the stator and rotor blade rows were separately 
calculated and the performance of the stage was also deter- 
mined. The rotor speeds varied from 9,500 to 18,600 r.p.m. 
and the pressure ratios varied from 1.75 to 3.25. The Mach 
numbers at the stator exit varied from 0.79 to 1.38. The 
results for this turbine are appraised and a procedure is 
demonstrated for smoothing loss coefficient data from the 
turbine rig . Test rig improvements reported include the 
design and construction of a new flow nozzle and the revision 
of the data reduction programs to access a Hewlett-Packard 
Model 9867B Mass Memory unit. 
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I. 



INTRODUCTION 



The transonic turbine test rig in the Turbopropulsion 
Laboratory at the Naval Postgraduate School was designed 
to determine the effect of blading design on turbine effi- 
ciency and to allow the separate determination of rotor and 
stator losses. 

Preceding the work of Solms [Ref. 5] in which several 
different turbine geometries were tested, the separation of 
the losses in the blade rows was not achieved reliably. 
Progressive improvements in the hardware and instrumentation 
of the turbine test rig following each series of tests, and 
finally the simplification of the data reduction program, 
resulted in the successful separation of losses for particular 
configurations reported in Reference 5. However, anomalies 
remained to be explained, and the values obtained for the 
stator and rotor loss coefficients were scattered. 

The questions raised by the preceding test program mainly 
involved the results of the turbine configuration designated 
in Reference 5 as Turbine C. Turbine C had converging- 
diverging stator passages in an axial entry, single impulse 
stage which was designed to operate at a pressure ratio of 
3.71. Results obtained for this geometry discharging to 
atmosphere did not agree with results obtained at reduced 
pressure levels when plotted in terms of similarity varia- 
bles ("referred" quantities) . Furthermore, the values 
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obtained for flow rate were of questionable accuracy since 
the existing flow nozzle was too large. Also the correction 
for a labyrinth leak rate became a more significant fraction 
of the turbine flow rate. 

The goal of the present work was to determine the blade 
row and stage performances of Turbine C, while resolving the 
anomalies reported in Reference 5. This would allow the use 
of the test rig to investigate effects of geometrical changes 
(stator-rotor separation and tip clearance for example) on 
the losses. The first step was to design and construct a new 
flow nozzle for the range of flow rates required by Turbine 
C. The flow nozzle was carefully calibrated against three 
different nozzles operating choked. The leakage rate for 
the labyrinth seal was also measured and a semi-analytic 
representation of the leakage rate as a function of pressure 
ratio was obtained which should apply at different tempera- 
tures. These developments are reported in Appendix A. The 
existing data reduction program [Ref. 5] was revised to 
access a Hewlett-Packard Model 9867B Mass Memory unit. The 
revised program is described fully in Appendix B. 

In the results of the test program given here in Section 
III, the effect of pressure level reported in Reference 5 
was not observed. The scatter in the loss coefficient however 
was not removed by the improvements in the accuracy of the 
measurements. In the discussion in Section IV, a method 
for smoothing the loss coefficient data is demonstrated 
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which will then allow a useful study of the effect of axial 
and tip clearances to be carried out. 
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II. TURBINE TEST RIG INSTALLATION 



A. DESCRIPTION 

The test installation consists of three major components: 
an Allis Chalmers twelve stage axial flow compressor, an 
exhauster assembly, and the turbine test rig (TTR) itself. 

The compressor is the source of driving air for the TTR 
and for the exhauster assembly. Fig. A-1 shows the piping 
arrangement. Turbine air passes through the first settling 
tank into an eight- inch pipe containing a flow nozzle, into 
the second settling tank and into the turbine. 

Fig. 1 shows the plenum, the floating stator assembly, 
the rotor, and the dynamometer [Ref. 5] . Pressure 'ratios of 
6:1 can be achieved when the system is hooded. The hood 
was needed to achieve near design pressure ratios in the 
tests reported here. Fig. 2 shows the turbine blading 
arrangement of the stator and rotor. Ref. 8 contains detailed 
description of the test rig hardware. 

The floating stator assembly depicted in Fig. 3 pennits 
measurements of the stator torque while axial and rotational 
movements are detected by calibrated force transducers that 
are heat insensitive. These measurements allow the deter- 
mination of the axial and tangential velocity components 
at the stator exit [Appendix A, Ref. 5] . 

In this report one configuration designated Turbine C 
was tested, the geometry of which is shown in Fig. 4; Table 
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1 describes the geometry quantitatively. The stator blade 
profile is shown in Fig. 5. The blades generate a converging- 
diverging nozzle shape. Pressure measurements are taken at 
the locations shown in Fig. 5. The pressxires necessary to 
the analysis of stator forces (Appendix B) are taken at the 
locations shown in Fig. 6 

B. TEST MEASUREMENTS AND ACCURACY 

MASS FLOW RATES : Appendix A gives a detailed analysis 

of the turbine flow rate and the labyrinth seal leak rate, 
and the results of calibration measurements. The results 
for the discharge coefficient are given in Fig. A-6. Fig. 

A-10 shows the Kinetic Energy Factor (Appendix A) obtained 
in labyrinth seal leak rate calibration and its respective 
polynomial. The equations used in calculating the flow rate 
from measurements are summarized in Table A-1 and Table A-2. 

FORCES, TORQUES, TEMPERATURES AND PRESSURE ; References 
8 and 9 give calibration procedures for the TTR. Identical 
procedures were employed here. Table II of Ref. 5 gives the 
expected accuracies of these measurements. 

C. TESTING AND DATA REDUCTION 

The TTR data collection system is described in Ref. 11. 
Appendix D in Ref. 5 gives a detailed explanation of the 
turbine test procedures. The data reduction program developed 
in Ref. 5 was revised to use the mass memory system. The 
revised program is described in Appendix B. 
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Seven tests were conducted of Turbine C at pressure 
ratios from 1.75 to 3.25. The parameters and conditions 
of the tests are summarized in Table II. 
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III. RESULTS 



A. INTRODUCTION 

The reduced data from seven tests of Turbine C are 
listed in Table V to XI. Table IV gives a listing of the 
symbols and meanings of the Table headings used in Table V 
to XI. Plots of the reduced data are presented in Figures 
7 to 16. Table III gives an explanation of symbols used in 
these figures. 

Tests 1 and 2 were conducted with the hood off. Tests 
3, 4, 5, 6 and 7 were made with the hood installed. Tests 
5 and 6 were repeats of tests 3 and 4 respectively. The 
results are discussed in the following paragraphs. ' 

B. TORQUE 

Fig. 7a and 7b give the results of the referred torque 
versus referred speeds at a constant pressure ratio. A 
nearly linear dependency exists for all seven runs with 
consistency in performance for the 1.75 and 2.25 hooded and 
unhooded pressure ratios. 

The results from the unhooded tests (Fig. 7a) agree with 
the results of the hooded tests (Fig. 7b) , at the same 
pressure ratio, to within 1 percent. 

C. EFFICIENCY (TOTAL-TO-STATIC) 

Figures 8a and 8b give the results of the total-to-static 
efficiency versus referred R.P.M. Figures 9a and 9b show 
plots of efficiency versus isentropic head. The maximum 
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efficiency occurred at similar values of the isentropic 
head for all pressure ratios. Ref. 17 relates efficiency 
to isentropic head. The peak efficiency was achieved at an 
isentropic head coefficient between 4.0 and 4.5 for both the 
hooded and unhooded tests. 

The efficiency is seen to increase as the pressure 
ratio is increased, approaching 80 percent at the pressure 
ratio of 3.25. The efficiency would be expected to increase 
to the design pressure ratio of 3.71, at which the stator 
nozzles should be just correctly expanded. 

D. THEORETICAL DEGREE OF REACTION 

Figures 10a and 10b give the results of theoretical 
degree of reaction versus isentropic head coefficient at 
different pressure ratios. Ref. 5 and Ref. 17 give the 
equation for the theoretical degree of reaction. For all 
runs the theoretical degree of reaction was negative at 
lower speeds. For a negative value, P 2 /P^q must be greater 
than values become positive only at high r.p.m. 

in both the hooded and unhooded configurations, passing 
through the design value of zero at a value of the isentropic 
head coefficient corresponding to the maximum in efficiency. 
Thus , negative values of the theoretical degree of reaction 
occur at speeds lower than the design speed at a given 
pressure ratio. Typically, there is seen to be also a loss 
of three to five percent in efficiency from the maximum 
value as the isentropic head coefficient is increased 100 
percent above the design value. 
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Results from hooded and unhooded tests are seen to 



agree with the exception of Run 3 in Figure 10b for which 
the values are low. Lack of agreement for the data of this 
test was found in all but the overall stage performance 
measurements . 

E. EFFECTIVE DEGREE OF REACTION 

Figures 11a and 11b give the results for the effective 
degree of reaction versus isentropic head coefficient. At 
all pressure ratios, hooded and unhooded, the effective 
degree of reaction was negative at lower speeds. Ref. 17 
states that the effective degree of reaction is a means for 
judging whether the flow is being accelerated or decelerated 
in the rotating cascade. The negative degree of redaction 
in Fig. 11a and 11b show that the flow is decelerated at 
lower r.p.m. and becomes accelerated as r.p.m. is increased 
(Isentropic head is decreased) . In no case should the design 
value of the effective degree of reaction be less than zero. 
As stated above, operation at speeds less than design speed 
results in a loss in stage efficiency, partially due to the 
undesirable deceleration in the rotor. 

The results for Run 4 (at a pressure ratio of 3.25) in 
Figure 11b reflect the marked improvement in efficiency at 
pressure ratios approaching design. 

F. STATOR LOSS COEFFICIENT 

The results for the stator loss coefficient are shown 
in Fig. 12a and Fig. 12b. While there is some scatter in 
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the data, the results from hooded and unhooded tests are 
seen to be in reasonable agreement. There is a slight 
increase in the stator loss as speed is increased at any 
pressure ratio. The stator loss for the highest pressure 
ratio (3.25) is significantly lower (13-17 percent) than for 
the lower pressure ratios (18-25 percent). This reflects 
the poor performance of a converging-diverging nozzle at 
pressure ratios less than design. 

G. ROTOR LOSS COEFFICIENT 

The results for rotor loss coefficient are shown in 
Figure 13a and Figure 13b. For completeness, the results of 
Run 3 have been included where the rotor loss coefficients 
were evaluated to be positive. Negative loss coefficients 
were obtained at higher speeds in Run 3, indicating an 
error in at least one element of the primary data. 

The rotor loss coefficient is the performance pareimeter 
most sensitive to small errors (and possibly unsteadiness) 
in the measurements. The sensitivity arises because it 
involves most of the parameters (see Ref. 5) previously 
calculated in the form of a ratio of small differences. In 
view of this, if the results of Run 3 in Fig. 13b are excluded 
from consideration, the rotor loss coefficient for the tests 
are reasonably consistent. The results for hooded and 
unhooded tests agree to within 20 percent at corresponding 
pressure ratios. There is also a definite decrease in the 
rotor loss coefficient as pressure ratio is increased toward 
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the design value, and a consistent decrease as the speed 
is increased at a fixed pressure ratio. The simultaneous 
reduction observed in the rotor and stator loss coefficients 
at higher pressure ratios is consistent with the increase 
measured in the stage efficiency (Fig. 8 and Fig. 9) . 

H. FLOW ANGLES 

An example of the variation in the flow angles at the 
stator and rotor exits as a fiinction of the isentropic 
head coefficient is shown in Fig. 14. As was found in all 
cases tested in Ref. 5, the variations in stator exit flow 
angle and rotor exit relative flow angle are very small. 

Also, the total variation of the relative flow angle into 
the rotor is only about 6 degrees for a 50 percent change 
in r.p.m. The results in Fig. 14 are qualitatively repre- 
sentative of the results for each of the seven tests. The 
exact variations in the measured angles can be seen in 
Tables V to XI. 

The stator exit flow angle is shown in Fig. 15a and 
Fig. 15b as a function of isentropic head coefficient for 
the seven tests. It can be seen that the measurements are 
in the range of 75 to 77 degrees, with the lower value 
occurring at the highest pressure ratio (3.25). It can be 
seen in Fig. 5 that the suction surface of the stator blade 
is designed to exit at 75 degrees, and the pressure surface 
at a smaller angle. Thus the measured stator exit angle 
appears to be slightly larger than would be expected; however. 



21 



the measured values when the nozzle passages are nearly 
correctly expanded approach the geometrical angle. 

The large values of stator exit flow angle at lower 
pressure ratios, and the higher values of stator loss 
coefficient are probably the results of the presence of flow 
separation. The effect of flow separation would be reflected 
in two ways. Firstly, the flow angle of the separated flow 
could depart significantly from the trailing-edge angle. 
Secondly, the existence of a region of separated flow could 
seriously affect the averaging process which is implicit 
in the data reduction. 

I . FLOW RATE 

The referred flow rate is shown as a function of the 
isentropic head coefficient for the seven tests in Fig. 

16a and Fig. 16b. The referred flow rate was constant at 
1.02 Ibs/sec to within the accuracy of the measurements 
for all speeds and all pressure ratios when operating 
hooded (Fig. 16b) . During unhooded operation (Fig. 16a) , 
the results gave a constant value of 1.01 Ibs/sec. 
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IV. DISCUSSION 



A. RESOLUTION OF PREVIOUS ANOMALIES 

In the results presented above, the difference in the 
turbine performance operating hooded compared with unhooded, 
is seen to be small. This is in contrast to the results pre- 
sented in Reference 5. Three factors contributed to the 
resolution of this anomaly; the accurate measurement of the 
turbine flow rate, installation of temperature insensitive 
load cells within the hood to measure stator loads, and 
corrections that were made to the data reduction program. 

There remains to be explained a difference of 1% in the 
flow rate between hooded and unhooded operation. More data 
are needed from unhooded tests to confirm that this differ- 
ence is repeatable. Because of the sensitivity of the perfor- 
mance parameters to the flow rate, the residual disagreement 
between hooded and unhooded results might be resolved if 
the difference in flow rates were explained. 

B. FLOW. RATE 

The use of a choked nozzle to calibrate a large flow 
nozzle was successful (Appendix A) . With the blockage 
factor correction applied, the results for different sizes 
of the choked nozzle agreed in the range for which the flow 
rates overlapped. 

The labyrinth leak rate calibration measurements were 
well represented by the analysis given in Reference 7. The 
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value obtained for the Kinetic Energy factor was higher than 
values given in Reference 5; however, the geometry of the 
turbine testing labyrinth was not known exactly. There was 
some eccentricity in the alignment of the stator section 
within the gland, so that the clearance was nonuniform to 
an unknown degree . 

The flow rate through the turbine, when expressed as 
a referred quantity, was measured to be a constant over the 
range of pressure ratios tested. (However, there was scatter 
in the measurements that resulted from lonsteadiness in the 
measurements of nozzle pressure drop.) Since the stator 
was choked, the constant value was expected. The value 
of the constant was consistent with the flow rate computed 
for the stator using a calculated value for the stator 
throat area. A value of the blockage factor within 1% of 
unity was required to be applied to the geometrical stator 
area to obtain total agreement. Since the exact area of the 
throat where choking occurs in the stator is not known to 
this accuracy, it is concluded that flow nozzle and stator 
throat measurements of the flow rate are in agreement. 

C. LOSS COEFFICIENTS 

The scatter observed in the results for the loss coeffi- 
cients makes a comparison between predicted and measured 
values inconclusive. There are trends in the data which are 
discernible, the rotor losses are consistently higher than 
the stator losses, for example. However, the scatter does 
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not allow the effect of parameter changes on the separate 

loss coefficients to be determined. The scatter is the 

result of the extreme sensitivity of the loss coefficient to 

the separate measurements on which they depend. For 

example, the stator loss coefficient (5 ) is given by 

s 






s 




( 1 ) 



where is the non-dimensional velocity and p^^ is the 
calculated pressure ratio station 1. Differ- 

entiation of Eq. (1) gives, if is constant. 
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Taking the sonic condition of P^^ = 0.532, with y = 1.4 and 
5 = 0.2, the quantity in brackets has the value 5.78. Thus 

o 

a 1 % error in calculating P^^ produces - 6 % error in the stator 
loss coefficient. Since P^ is derived (as is X^^) from 
measurements scanned over a period of more than a minute, 
with some variation in operating point (as well as an 
observable fluctuation in the flow nozzle pressure drop) 
occurring during the data scan, the scatter in the loss 
coefficient data is understandable. 
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Before the loss coefficient can be used as a measure 



of performance, a method of smoothing the data while 

retaining accuracy is needed. Two attempts were made 

here. First, an attempt was made to smooth the calculated 

values of the interstage pressure (p^) , by comparing p^^ 

to the pressures measured at the hub and tip. Figure 17 

shows the variation p^ and the measured pressure with 

isentropic head coefficient for one particular test. 

Clearly, the calculated average pressure at station 1 must 

be smoothly behaved if the hub and the tip pressures are 

smoothly behaved. An attempt to derive a smoothing function 

for the interstage pressure is shown in Fig. 18, where the 

data from several tests at the same pressure ratio are 

plotted in a dimensionless form. The figure clearly shows 

the degree of scatter in the values of P^. On the basis of 

these data above, however, the use of a single polynomial 

for the quantity (P,-P, i_)/(P^- -P. . ) as a function of 

^ -' I hub tip hub 

isentropic head coefficient cannot be justified. When a 
polynomial was derived from the data of a single test, and 
the loss coefficients for the stator and rotor were recomputed, 
the results shown in Fig. (19a) and (19b) were obtained. 

The second attempt to smooth the loss coefficient 
data involved effectively the elimination of scatter in the 
input data from which the interstage pressure and velocity 
were calculated. The referred flow rate was determined to 
be constant when the stator was choked. By setting the 
referred flow rate equal to the constant value (1.02 or 1.01, 
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hooded and unhooded, respectively) and continuing the data 
reduction as before, the results shown in Fig. 20a and 20b 
for the loss coefficient were obtained. A considerable 
reduction in scatter is evident, with the major benefit in 
the rotor loss coefficient. Further data must be examined 
to fully appraise the benefit of this smoothing technique. 

A third technique, so far not attempted, is to calcu- 
late in the data reduction before calculating the axial 
velocity component at station 1. The values of could 
be smoothed using a polynomial curve fit and the smoothed 
values of P^ used to calculate the velocity. An application 
of all three techniques should produce loss coefficient 
distributions which are smooth and yet accurate. 
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V. CONCLUSIONS AND RECOMMENDATIONS 



The following conclusions are made for turbine configura- 
tion C. 

1. The loss coefficients are extremely sensitive to 
input data, especially the rotor loss coefficient. Repeata- 
bility in calculating the loss coefficient directly from 
measurements was poor. 

2. Methods were demonstrated to smooth the loss coef- 
ficient data without sacrificing, and possibly improving, 
the accuracy. 

3. For all pressure ratios there was a velocity 
deceleration through the rotor and a pressure rise 'across 
the rotor at lower r.p.m., indicating perhaps a shock wave 
at or in the rotor. 

The following recommendations are made: 

1. Apply the three data smoothing techniques described 

above to all data from Turbine C in sequence: input constant 

referred flow rate, calculate interstage pressure and smooth 
using measured hub and tip pressures, then calculate 
interstage velocity. 

2. Determine a more accurate input for K^, the ratio 

of the momentum average velocity to the mass-average velocity 
and examine the effect on the calculated stator exit flow 
angle. 
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3. Program the quadratic curves for the theoretical 
loss coefficient in references 16 and 18 in order to reduce 
the input error. 

Finally, it is concluded that following the recommenda- 
tions listed here, the effect of parameter variations on 
the blade row losses could be examined satisfactorily using 
the turbine test rig. 
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TABLE I 



TURBINE GEOMETRy 





TURBINE 


STATOR 


ROTOR 


C 


1 


1 


DESCRIPTION 


A* 2 

(in) 


I NUMBER 
OF 

! BLADES 


“XIN) 


h/ . >, 

(in.) 


STATOR 1 


CONVERGING 

DIVERGING 

NOZZLES 


2.9058 


51 


4.184 


0.5775 


ROTOR 1 


CIRCULAR 

ARC 


7.119 


60 


R„1-4.195 

R^2=4.250 


h^=0.752 

h2=0.8475 
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TABLE II TEST PARAMETERS FOR TURBINE C 
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TABLE IV 



DEFINITION OF COLUMN HEADINGS IN TABLES V. TO XI 



SYMBOL 


DEFINITION 


VI 


VELOCITY AT STATOR EXIT PLANE (FT/SEC) 


V2 


VELOCITY AT ROTOR EXIT PLANE (FT /SEC) 


VAl 


AXIAL VELOCITY AT STATOR EXIT PLANE (FT/SEC) 


VA2 


AXIAL VELOCITY AT ROTOR EXIT PLANE (FT /SEC) 


VUI 


TANGENTIAL VELOCITY AT STATOR EXIT PLANE (FT/SEC) 


VU2 


TANGENTIAL VELOCITY AT ROTOR EXIT PLANE (FT/SEC) 


MI 


MACH NUMBER AT STATOR EXIT PLANE 


MAI 


AXIAL MACH NUMBER AT STATOR EXIT PLANE 


M2 


MACH NUMBER AT ROTOR EXIT PLANE 


MA2 


AXIAL MACH NUMBER AT ROTOR EXIT PLANE 


A1 


FLOW ANGLE AT STATOR EXIT PLANE (DEGREES) 


A2 


FLOW ANGLE AT THE ROTOR EXIT PLANE (DEGREES) 


BI 


RELATIVE FLOW ANGLE AT ROTOR INLET PLANE (DEGREES) 


B2 


RELATIVE FLOW ANGLE AT ROTOR EXIT PLANE (DEGREES) 


ZS 


STATOR LOSS COEFFICIENT 


ZSTH 


THEORETICAL STATOR LOSS COEFFICIENT 


ZR 


ROTOR LOSS COEFFICIENT 


ZRTH 


THEORETICAL ROTOR LOSS COEFFICIENT 


ZR* 


ROTOR CARRY OVER LOSS COEFFICIENT 


ZI 


ROTOR INCIDENCE LOSS COEFFICIENT 


P.R. 


PRESSURE RATIO P /P 
to 2 


K-IS 


ISENTROPIC HEAD COEFFICIENT 


NREF 


REFERRED RPM (R.P.M.) 
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WREF 


REFERRED FLOW RATE ( LBM/SEC ) 


ETA 


TOTAL TO STATIC EFFICIENCY 


HP REF 


REFERRED HORSEPOWER ( H.P. ) 


RTH 


THEORETICAL DEGREE OF REACTION 


REFF 


EFFECTIVE DEGREE OF REACTION 


PI 


PRESSURE RATIOS CORRESPONDING TO TAP LOCATIONS 


P2 


IN FIGURE 5. PRESSURES ARE REFERRED TO STATOR 


P3 


INLET TOTAL PRESSURE 


P4 




P5 




P6 




P7 




P8 




P9 


V 


KBLOCK 


STATOR BOCKKAGE FACTOR 


STPR 


PRESSURE RATIO ACROSS THE STATOR P /P, 

to 1 


RPM 


ROTOR SPEED ( R.P.M. ) 
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FIGURE li TRANSON'IC TURBINE TEST RiO 
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FIGURE 4 TURBINE C 
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FIGURE 6 TURBINE TEST RIG GEOMETRY FOR 
TURBINE GONE IGUR.^T ION C 
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FIGURE 7a REF. ROTOR TORQUE VS REF. SPEED 
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FIGURE 7b REF. ROTOR TORQUE VS REF. SPEED 
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FIGURE 8a EFFICIENCY VS REF. SPEED 
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FIGURE 8b EFFICIENCY VS REF. SPEED 
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FIGURE 9a EFFICIENCY VS . ISENTROPIC HEAD COEFFICIENT 
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FIGURE 9b EFFICIENCY VS ISENTROPIC HEAD COEFFICIENT 
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FIGURE 10a TH. DEGREE OF REACTION VS HEAD COEFFICIENT 
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FIGURE lOb TH. DEGREE OF RACTION VS HEAD COEFFICIENT 



i 1 F 



^ 1 1 1 ^ 



+ 



+ 



ui 

m 



UJ 

Cx 

a 

a 



E 

B 

LJ 

VC 

Er^rvjLn 
• • 

Q.— fMO; 

cu 

— TJU. 

a 



tranu 

»»>• 

C3XLU 



□3 



m 



X 

CO 

X 

X 

D 

X 



aif 



X 

cn 



X 

DO 



X 

cn 



r- 



1/1 

La 



H 

:3 

u 

M 

a 

M 

Pm 

Pm 



Ln 

1/1 



Ln 



Ln 

m 







< — ± — I — ^ — I — i — i 

nj ca rvj 

• • • 

s s s 

I 

ETFCaiVE DEG. OF 



A I A 

X to 

• • 

es es 

1 i 



REflCTIDN 



m 

• 

cs 

i 



a 

B 



111 

IM* 



I 



61 



15ENTRGPIC HERD CPEFFKIENT 
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15ENTRDPIC HERD CDEFFICIENT 

FIGURE 13a ROTOR LOSS VS. ISEOTROPIC HEAD COEFFICIEHT 
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15ENTRDPIC HERD CDEFFICIENT/HDDDED 

FIGURE 13b ROTOR LOSS VS. ISENTROPIC HEAD COEFFICIENT 






] 



1 00 . 001 



H — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I- 



alb X 

^ X 

m X 
m X 

m X 
^ X 
♦ X 
m X 



in 

rvj 

ni 

a 

• 

VC 

CL. 



CC 

cx 

Ui 

cx 

in 

— rvi — pj {2 

CC CC £ 

zc n: CC CC 
GL d. I— I— u. 
_J ^ LiJ LU CC 
£E CC tn tfl 

n n ti It ^ 

t4. 

4% X a b 



m 

m 

m 



m 



m 



cn 



1/1 

m 



1/1 



t/i 

in 



LlJ 

iCs 

a 

a 



LiJ 

C 



UJ 

t/1 C3 

• 

1/1 

CC 



1/1 

3" 



CU 

a 

cx 



LU 

in 



1/1 

n 



1/1 













1 ■« cb ' 






a pi 


S 


ca 


63 


s 


s 


ca B 


s 


B 


S 

m 


63 

to 


63 

X 


63 

nj 


s 


S 63 

rj X 


63 

U3 


63 

m 



TLOM angle: degrees 



67 



FIGURE 14 FLOW ANGLE VS. ISENTROPIC HEAD COEFFICIENT 
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FIGURE 15a FLOW ANGLE ( ALPHA 1 ) VS. ISENTROPIC HEAD COEFFICIENT 
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FIGURE 16a REF. FLOW RATE VS. ISENTROPIC HEAD COEFFICIENT 
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I5ENTRDPIC HW COEFFICIENT 

FIGURE I6b REF. FLOW RATE VS. ISENTROPIC HEAD COEFFICIENT 
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FIG, p calculated ani> measured interstage pressures 
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FIG. 19a EFFECT OF SMOOTHING INTERSTAGE PRESSURE ON STATOR LOSS COEFFICIENT 
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FIG. 19b EFFECT OF SMOOTHING INTERSTAGE PRESSURE ON ROTOR LOSS COEFFICIENT 






! ■ I - I— h i I I - r I I I I I I -I H ■ !■ I I 



t/i 

o 



XQl 



Lrt 



H* Ui. 

0 Zl 

a I— 

S §' 

1 s: 

Lrt Lrt 
W W: 

ni fU 

UK; uti 



ES 

z as. 

I I 
C 3 - x: 



XOK 



Lrt 

m 



lrt 

m 

Lrt 



Lrt 



LlI 



hi 



Ld 

CZL 



LxJ 

in 



ut 

m 



m n ni rvi — — 



t/r 
ru 



5THT0R LD55 CDEFFICIENT 



76 



FIG, 20f EFFECT OF ASSUMING CONSTANT REFERRED FLOW RATE ON STATOR LOSS COEFFICIENT 
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APPENDIX A 



FLOW RATE DETERMINATION 



A-1 INTRODUCTION 

A precise measurement of the flow rate is necessary 
if the performance of an individual blade row is to be 
determined acctirately from measurements made with the turbine 
test rig. Anomalies that occurred in the results from the 
first tests of the present supersonic turbine (Turbine C 
in Ref. 5) were possibly due to inaccuracies in the speci- 
fied flow rates. Since the flow rate for Turbine C was 
smaller than for the turbines previously tested in the rig, 
small pressure differentials were measured with the existing 
nozzle. Also, the higher pressure ratios at which Turbine 
C was operated resulted in increased leakage through the 
labyrinth seals, the geometry of which might have change 
slightly on reassembly. Therefore, before new tests of 
Turbine C were begun, a new flow nozzle was provided and 
careful measurements were made of the labyrinth leak rate. 

The flow nozzle design and calibration are described in 
Section A-2. The labyrinth leak rate determination is given 
in Section A- 3. 

A- 2 MEASUREMENTS OF THE TOTAL FLOW RATE 

A-2.1 FLOW NOZZLE DESIGN 

The flow nozzle is positioned in a pipe through which 
air is supplied to the turbine (Fig. A-1) . The nozzle 
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generates a differential pressure to which the flow rate 
is related. 

For a one-dimensional steady flow, the mass flow rate, 
W, is given by 



W = p AV 



(A-1) 



where p is the density, V is the velocity, and A is the 
cross sectional area. 

Using the perfect gas equation of state 



P 




(A-2) 



and Bernoulli equation 



P 



t 




(A- 3) 



Eq. (A-1) for low velocities becomes 



W = A 






(A-4) 



The symbols P and T are pressure and temperature, respec- 
tively, and the subscript t denotes stagnation conditions. 
This equation, with explicit corrections for compressibility 
and thermal expansion, is used to calculate the flow rate 
through a flow nozzle from measurements of temperature and 
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pressures. Withe the addition of the correction factors and 
with conversion factors added to account for the units of 
measurements [Ref. 1] , Eq. (A- 4) becomes 



W = 
N 



0.16384 

/2.036 



N 






ij 



(Ibs/sec.) (A-5) 



where 



N 

Y, 



N 

Ah 



0.16384 



Diameter of nozzle (inches) 
Thermal expansion coefficient 
Compressibility coefficient 
Pressure at the nozzle (psia) 
Water differential (in. H 2 O) 
Temper a tur e ( ° R) 

Conversion factor 



and lastly, is the discharge coefficient. 

The discharge coefficient accounts for the fact that the 
flow is not one-dimensional, and that the pressures measured 
are not and (P^ - P) . Generally, the discharge coefficient 
has a value above 0.95 for a standard ASME nozzle at high 
Reynolds numbers. Establishing the relationship between 
Reynolds mamber and the discharge coefficient is what con- 
stitutes the calibration. 
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The flow rates expected in the supersonic turbines were 
1.3 to 3.5 Ibm/sec. A flow nozzle with a diameter of 3.25 
inches was chosen (by applying Eq. (A-5) with unity for the 
coefficients) to give sutiable water column measurements of 
approximately 9.8 to 58.1 inches. The nozzle design followed 
the criteria given in Ref. 1 for "Low Beta series" designs. 
The geometry of the flow nozzle is shown in Fig. A-2 and 
views are given in Fig. A-3. The locations of the pressure 
taps are shown in Fig. A-2. Four throat taps are located 
1.5 nozzle diameters from the nozzle face and are manifolded 
together. The upstream pressure is measured using taps in 
the upstream flange. 

A-2. 2 FLOW NOZZLE CALIBRATION 

The discharge coefficient is given in terms of measure- 
ments, if the flow rate is determined by independent means. 
One of the most accurate methods of determining the flow 
rate of a gas is by choking the flow at a known area. This 
is because the mass flux is nearly xaniform where choking 
occurs, and because the boundary layer can be kept thin in 
a rapid smooth contraction from a relatively large pipe. 
Hence, in order to determine the discharge coefficient of the 
new flow nozzle as a function of Reynolds number, the 
arrangement shown in Fig. A- 4 was used. 

A-2. 2.1 Method 

Flow nozzle calibration runs were made at con- 
trolled supply pressures from 15 psia to 45 psia. Three 
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different orifice plates with diameters of 1.6, 2.065, and 
2.24 inches were mounted in ttirn to provide choked flows 
over the desired range of Reynolds numbers. The differen- 
tial pressiire at the flow nozzle was varied in increments 
of 2 to 4 inches of water. The range of flow rates over 
which the orifice plates were choked at the supply pressure 
was from 1.33 to 3.55 Ibm/sec. The maxim^lm pressure ratio 
obtained was 2.92. 

A-2.2.2 Analysis 

The discharge coefficient was given by Eq. (A- 5) , 
and the Reynolds number was determined by 






48 Wjj e 



7T D y 
P 



(A- 6) 



where y is the viscosity and 6 is the ratio of the diameter 
of the nozzle to the diameter of the pipe. The flow rate, 
which was the same for the nozzle as for the choked plate, 
was determined by the following equation 



r ^t ^ S 

/RT^/g 



(A- 7) 



where 



= pressure at choked nozzle (psia) 

2 

A = area of choked nozzle (in ) 
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Kg = blockage factor 
R = gas constant {ft-lb^/lbm - °R) 

= Temperature at choked nozzle (®R) 

The function r is a known function of the specific heats and 
is given by 



The blockage factor Kg accounts for the boundary 
layer and for the axisyinmetric surface; Kg is given by 
[Ref. 3] 



K 



B 



1 - 



45* 

°*N 



(A-9) 



where the star on signifies sonic conditions, and 5* 

is the boundary layer displacement thickness. 

The boundary layer growth was estimated using 
the method of Falkner and Garner given in Ref. 6. The 
displacement thickness, 5*, is given by 



* 

5 




Ha 



n 

n+1 



Re 



1 

n+1 

1 



n 



' 1 ^ 



n+1 



(A-10) 
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where is the distance from the stagnation point to the 
choking station. H is the form factor, and a and n are 
empirical constants . i|; is a factor which depends on the 

pressure gradient along the surface. Re^^ is the Reynolds 
nvimber based on conditions at the throat and s^^. 

It is shown in Ref. 6 that the value of ip for 
a similar contraction was 0.38. With \p = 0.38, and with 
H = 1.4, a = 0.0076 and n = 6 for air, Eq. (A-10) becomes 



6 * = 



0.0932174 



Re 



0.1429 "1 



The Reynolds number is given by 



W St 

Re, = 12 ^ — 
1 



(A- 11) 



(A-12) 



where is the viscosity at the throat where the tempera- 
ture is 5/6 times the stagnation temperature. 

For each operating condition for each orifice 
plate, Re^^ was calculated from Eq. (A-12) using W calculated 
from Eq. (A-7) with K„ = 1.0. 6* was then calculated using 

Eq. (A-11) and K„ obtained from Eq. (A-9) . W was recalcu- 
lated from Eq. (A-7) . It was found that the blockage factor 
could be taken to be constant with pressure for each choked 
plate at the following values: 
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PLATE if 


THROAT DIAMETER 


s^ (ins.) 




1. 


1.6 


3.51 


0.9928 


2. 


2.065 


3.10 


0.9935 


3. 


2.24 


3.22 


0.9901 



The flow chart for the analysis of the calibra- 
tion test results is given in Fig. A-5 and the equations 
used are summarized in Table A-1. 

A-2.2.3 Results 

The nozzle calibration extended over a range 

5 5 

of Reynolds n\ambers from 0.081 x 10 to 0.213 x 10 . The 
results are given in Fig. A-6, where the data for three 
different orifices are shown for the range of pressures 
for which the flow was choked. A second order polynomial 
was found to represent the flow nozzle coefficient as a 
function of the Reynolds number. 



= 1.0272 



Re. 



0.1598 



N 



Re 



( — T 
10 ^ 



) + 0.3805 ( 



2 



10 



(A- 13) 



where Re^ is given by Eq. (A-6) , was used to reduce data 
from tests of the turbine test rig. 

A-2.3 APPLICATION 

The discharge coefficient is a function of Reynolds 
number. However, in application the flow rate is to be 
determined and therefore the Reynolds number is unknown. 

The determination of the flow rate from measurements requires 
the solution of three equations which have the general forms: 
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(A-14) 






(A-15) 



Re^j CC2 Wjj 



(A:=16) 



These equations represent Eq. (A-5) , Eq. (A-13) , and Eq. 

(A- 6) , respectively, where and U 2 are known in terms of 
the geometry and measurements for any operating point. 
Clearly, by substituting for Re^^ and using Eq. (A-15) , 
and Eq. (A-16) in Eq. (A-14) , a single function for the 
unknown W is obtained. 

An iterative procedure is used to solve the three 
equations for the unknown flow rate. The procedure is 
given in Appendix B. 

A- 3 MEASUREMENT OF THE LABYRINTH LEAK RATE 

A- 3 . 1 Labyrinth Leak Rate Calibration Test 

Since the leak rate of the plenum labyrinth cannot be 
measured during turbine operations, the relationship of the 
leakage flow rate to supply and hood pressures and supply 
temperature was determined in a separate calibration test. 
For this test, the main flow through the turbine was blocked 
by a plate and the labyrinth was allowed to leak into the 
hood as in normal operation. The flow was passed through a 
two- inch pipe (Fig. A-4b & Fig. A-1) , containing an ASME 
standard sharp-edged orifice. The usual supply pipe was 
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blocked. The flow was discharged through the labyrinth 
seal into the hood in which the pressure was controlled 
by the ejector. 

In determining the leakage rate using the ASME orifice, 
vena contracta taps were used. The differential pressure 
was measured in inches of water on a manometer which had 
0.1-inch graduations. The upstream pressures in the lab- 
yrinth and in the nozzle were measured in inches of mercury 
on U-tube manometers with atmospheric pressure as a refer- 
ence. The supply pressure was varied from 14 to 83 psia, 
while the differential pressure across the vena contracta 
taps varied from 2.75 to 37.3 inches of water. The maximum 
pressure ratio across the labyrinth was 5.6. 

A- 3 . 2 Calculation of the Leakage Flow Rate 

Since the flow orifice was in a standard installation 
according to the criteria given in Ref. 1, Table 5 of Ref. 1 
was used to obtain the discharge coefficient as a function 
of Reynolds number. In order to program the data reduction, 
the data of Table 5, for the appropriate values of Beta, 
were first approximated in the range of Reynolds numbers 
from 8,000 to 1,000,000 by a polynomial as shown in Fig. 

A-7. The discharge coefficient at each operating point was 
then given by 
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10 

K„ = 0.603698259 + 0.004688755 (i^) 
N RSj, 



(A-17) 



5 5 

- 1.250025 X lO"^ (4^)^ + 2.0320771882 x lO"^ ^ 



5 5 

- 1.56980769 x 10“^ (i^) + 4.5277098 x 10“^ (i^) ^ 



The flow rate was calculated at each point using Eq. 
(A-5) , using the expressions for the compressibility and 
thermal expansion coefficients given in Table A-2, and the 
discharge coefficient given by Eq. (A-17) . 

A-3.3 Analysis of the Results 
A- 3 . 3 . 1 Simple Method 

A non-dimensional flow rate, or flow function, 
for the flow through the labyrinth can be defined as 



4 > 



/(R/g) T^ 



(A-18) 



During the calibration tests , the stagnation 
pressure (P_) was varied in steps for constant values of the 

Li 

hood pressure • The stagnation temperature was approxi- 

mately constant. The flow ftinction was calculated for each 
operating point using Eq. (A-18) and the values plotted as 
a function of the pressure ratio results are 

given in Fig. A- 8. 

A second order polynomial was found to fit the 
data to an accuracy of ± 4% as shown in Fig. A- 8; 
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(A-19) 



$ = 0.2721 + 0.0365 (P, /P^ ) - 0.2357 (P^/P^ ) 

n L n L 



2 



While Eq. (A-19) represents the calibration tests 



results quite well, it is an empirical result for a limited 
variation in test parameters. Operation at different 
temperatures, for example, could give a different behavior. 



assumptions, the flow through a labyrinth with many teeth 
(N) can be represented by 



and Ke is the "Kinetic Energy Factor”. Here, the dependence 
of the leak rate on the upstream and downstream pressures 
(3) has been obtained analytically. The Kinetic Energy 
Factor is a measure of how much kinetic energy is recovered 
as stagnation pressure as the flow passes successive teeth. 
It is expected to be a function of Reynolds number, and 
hence would include the correct behavior with stagnation 
temperature. 



A-3.3.2 Use of the Kinetic Energy Factor 
It is shown in Ref. 7 that, under reasonable 



$ = 3 Ke 



(A-20) 



where 




(A-21) 
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The geometry of the labyrinth is given in Fig. 
A-9. Using the notation in Fig. A-9, the Reynolds mamber 
for the leakage flow is given by 






12 C 

Li 



(A-22) 



where 



A|^ = 2 IT D C (A-23) 

and li,. is the viscosity based on the labyrinth supply 
temperature, T,. . 

For the flow rates obtained in the calibration 
tests, $ was obtained from Eq. (A- 18) . 3 was calculated 

using Eq. (A-21) and Ke determined using Eq. (A-20) . The 
Reynolds number was obtained using Eq. (A-22) with Eq. (A-23) . 
The Kinetic Energy Factor was plotted as a function of the 
Reynolds number and the result is shown in Fig. A-10. A 
fourth order polynomial was found to fit the data to an 
accuracy of ± 2%: 

Ke = 0.601373889 + 4.36695141 x lo""^ (-^) 

10 



- 2.4530860 x lo"^ (-^) ^ + 6.5672233 x lO"^^ (-^) ^ 

10 10 

- 6.086573 X lO"^^ (^ 1 )"^ (A-24) 

10 
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It is noted that the use of Eq. (A-24) and Eq. 
(A-21) gives a factor of two improvement in accuracy over 
the use of the purely empirical result given as Eq. (A- 19) . 
Also, the dependence of the Kinetic Energy Factor on the 
supply temperature is at least partially included in the 
representation as a function of Reynolds number. 

The flow chart for the calibration data reduction 
is given in Fig. A-11 and a siammary of the equations is 
given in Table A-2. 

A-3.2 APPLICATION 

In order to use Eq. (A-24) to calculate the leak rate 
at a given operating point, an iterative procedure is 
necessary since the Reynolds number is again a function of 
the flow rate. The procedure is similar to that followed 
in calculating the total flow rate (Section A-2. 3) and is 
described in .^Vppendix B. 
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TABLE A- 


•1 




SUMMARY OF FORMULAS FOR 


NOZZLE CALIBRATION DATA REDUCTION PROGRAM 


FLOW 






0.16384 






• 

T.T 






r 


>/ \ \j 


I 




‘^N 




/ 2.036 ■ 




RATE 




> 










THEMAL 

EXPANSION 

COEFFICIENT 






1+0. 


00193 


T„ - 528 
N 

100 




COMPRESSIBILITY 

COEFFICIENT 


Y 

1 


= 


1 - ( 


D. 05246 


^ 0.41 + 0.35g^ 


Ah 

P 

N 


SETA 


0 




E 


N 


D = 3.25 ins. 

N 


, D = 7.975 ins 
P 




E 


•p 


VISCOSITY 


M 




1.153 X 10 X 0.06333 r, 

ri 










198.72 / T^ + 1 


• 


REYNOLDS 


73 






48 ^ ^ 




NUMBER; 


®N 






n D 

P 






FLOW 

FUNCTION; 


r 




k 


m 

2 


Y+v 

y-1 




Y+ 1 


BLOCKAGE 


K 

B 




1 « 








FACTOR; 












BOUNDARY 

LAYER 

THICKNESS 


X* 




0.21364 


(0.38 ) s 




0 








1 


REYNOLDS 

NUMBER;* 






1 ? - 


W 






RSi 












TI33i7 

RATE 


Us 




Pc A K B 


r 




CHOKED 

NOZZLE 


w 








I 








J RVg 
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TABLE A- 2 

SUMylARY OF FORMULAS FOR LABYRINTH CALIBRATION REDUCTION PROGRAM 



FLOW 


. Trt • 

y = 2.036 w,y 


FUNCTION: 


Pl 


FLOW 


0.16384 - rP Au' 

Wt = f D A„ Y, / — 


RATE: 




BETA: 


/ 1 - (!h)" 

^ N = 10 

V N - LS- (P^/ P^) 


REYNOLDS 


12 W, C 


NUMBERS: 


11 


AREA 

LAB: 


= 2*ir D C 


KINETIC 

ENERGY 

FACTOR: 


I 

Ke = £. 

(9 


THERMAL 

EXPANSION 

COEFFICIENT: 


Tn - 528] 

A = 1 + 0.00193 

" ■ 100 

• J 



COMPRESSIBILITY ^ h 



COEFFICIENT: 


^ ^ , 

^N 


BETA; 


Q 

P " * » - 3.25 ins., D„ - 7.975 ins. 

D P 

P 


VISCOSITY: 


1.153 X 10 ^ X 0.06333/t^ 
198.72 / Tj, + 1 
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FIGURE A-1 TEST RIG PIPING INSTALLATION 
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FIG. A- 3 VIEWS OF THE FLOW NOZZLE 
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FIG. A-4a PIPING ARELINGEMENT FOR LABYRINTH LEAKAGE CALIBRATION TESTS 
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FIG. A-4b PIPING ARRANGEMENT FOR FLOW NOZZLE CALIBRATION TESTS 
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FIG. A-5 FLOW NOZZLE CALIBRATION DATA REDUCTION FLOW CHART 
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FIG. A- 6 FLOW NOZZLE CALIBRATION 
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FIG. A- 7 ASME DISCHARGE COEFFICIENT 
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FIG. A- 8 PLENUM LABYRINTH SEAL LEAK RATE 
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FIG. A-9 PLENUM LABYRINTH 
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FIG. A- 10 LABYRINTH SEAL KINETIC ENERGY FACTOR 




FIG. A- 11 LABYRINTH CALIBRATION DATA REDUCTION FLOW CHART 
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FIG, A- 11 LABYRINTH CALIBRATION DATA REDUCTION FLOW CHART 
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APPENDIX B 



TURBINE TEST RIG (TTR) DATA REDUCTION AND PROCESSING 
B-1 INTRODUCTION; 

This appendix docioments the current data reduction 
program which was revised and updated from the program used 
in Ref. 5. The present program uses the "Mass Memory" disc 
storage unit attached to the Hewlett-Packard Model 9830A 
Calculator. The equations developed in Ref. 5 apply to 
the present program. 

A description of the programs is given in Section 
B-2 and the variables are listed and defined. Procedures 
for running the program are given in Section B-3. 

B-2 DESCRIPTION OF PROGRAMS; 

The reduction program is divided into eight separate 
sub-programs. Fig. B-1 shows the contents of each sub- 
program. Table B-1 gives the channel and port nxxmbers for 
each of the test measurements and the respective matrix 
elements assigned in the data reduction program. Table B-2 
defines the variables and equations (Ref. 5 & 16) used in 
the program. Table B-3 is a record or score sheet of the 
variables used in the program. 

B-3 PROCEDURE FOR DATA REDUCTION PROGRAM USING MASS MEMORY 
1. CHECK; 

a. ) Mass Memory to UNLOAD Position 

b. ) ROM installed 

c. ) NO cassettes in transport 
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2 . 



Turn on HP 9830A Calculator 



3. Turn on HP9866A Printer 

4. Turn on HP 9868A I/O Expander 

5. Turn on HP9863A Tape Reader 

6. Turn on HP 9867B Mass Memory 

7. Turn on HP 11305A Controller 

8. Load platter number TTR when the DOOR UNLOCKED light is 
lit. 

9. Switch to LOAD position; wait 30 seconds until DRIVE 
READY light is lit. 

10. Key in; GET ”TTR 1;” Wait for completion 

11. Press "RUN EXECUTE." 

12. The display will read "ENTER NEXT RECORD # ON DATA FILE." 
Input your next record number for raw data storage and 
press "EXECUTE." 

13. The display will read "TAPE: 1ST HOLE7-ON START? ;-CONT. " 

By pressing "CONTINUE, EXECUTE" the paper tape will be 
read. 

14. After the paper tape is read, the display will read 
"CORRECTIONS TO DATA? YES=1, NO=0;" This option permits 
you to correct any erroneous corrections. PRESS "1 
EXECUTE" if you have corrections. Press "0 EXECUTE" 

if no corrections; go to step 19. 

15. With Corrections to be made, the display will read 
"PRESS PRT ALL KEY FOR RECORD." This enables you to 
have a record of your data points. 
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16. The display will read "ENTER CORRECTION AS MATRIX 
ELEMENTS." When a scanivalve is detected to have given 
an erroneous value, the manometer board can be read 
and correct value recorded. See Table B-2 for the 
appropriate matrix identification. 

17. The display will read "ENTER CORRECT VALUE = ?;" 

Input the correct value and press "EXECUTE." 

18. The display will read "ANY MORE CORRECTIONS? YES=1, 

NO=0." Press "1 EXECUTE" if more corrections; go to 
step 17 . 

19. 1.5 inches of paper will emerge from the printer and 
the raw data will be printed. Afterwards the display 
will read "STORE DATE? ENTER YES=1, NO=0." If. raw 
data is to be stored in Mass Memory press "1 Execute." 

The raw data is stored. Below the tables will be 
printed the words "THE RAW DATA IS STORED IN RAW DAT 
FILE #? Along with the file n\unber that was input in 
step 12. If data is not to be stored press "0 EXECUTE" 
and the words "THIS DATA WAS NOT STORED" will be printed. 

20. The calculator will automatically CHAIN with "TTR 2." 

21. Temperatures of the nozzle, labyrinth and T^^ will be 
printed in degrees Rankine. The labyrinth seal flow 
rate, its Reynolds number and kinetic energy factor 
will be printed. Also, the nozzle discharge coefficient, 
nozzle flow rate, its Reynolds number, and flow rate 
will be printed. 
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22. The display will read "CALC NEW W-DOT? YES=1, NO=0." 

If you elect to calculate a new theoretical flow rate^ 
based on a stator blockage factor of 0.965, and employ 
this value in^ the remainder of the program, press "1 
EXECUTE." However, if you use the original flow rate, 
press "0 EXECUTE." Go to step 24. 

23. The election for a new calculation of flow rate was 
chosen; a new flow rate, the stator flow function and 
blockage factor will be printed out. Note that the 
flow function and blockage factor were employed to 
calculate the new flow rate. Go to step 25. 

24. The original flow rate is used. The stator flow function 
and blockage factor based on the flow rate and measure- 
ment are printed. 

25. The display will read "INPUT ALPHA 1 (DEG), OR 0." 

Input ing any number other than 0 will fix the stator 
exit flow angle in degrees. This method assvimes VU^^ 
is calculated from the stator torque and adjusts the 
value of to agree with the elected stator exit angle. 

By pressing "0 EXECUTE" the stator exit angle is not 
chosen. Go to step 27. 

26. By choosing to fix the stator exit flow angle, the 
value you selected will be printed. The value of K, 
which is the ratio of VAj^ determined by the momentum 
equation to VA^^ determined by continuity, will be 
printed. VAlues from 1.0 to 1.2 are within acceptable 
limits. Go to step 29. 
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27. The value of the stator exit as calculated from 
measurements will be used. The display will read "ENTER 
K." Values between 1.0 and 1.2 are within acceptable 
limits. Press " (K) EXECUTE." This value will be 
printed and the program will CHAIN with "TTR 3." 

28. TTR 3 will evaluate Control Volume B and will CHAIN 
with "TTR 4." 

29. The display will read "CALC THEOR. LOSSES ? YES = 1, 

NO = 0 . " This program option permits you to calculate 
both the theoretical rotor and stator loss coefficients. 
Copies of REF. 16 and 18 will be needed. By pressing 

"1 EXECUTE," you have elected to calculate the losses. 
Press "0 EXECUTE" if you do not desire to calculate 
the losses. Go to step 37. 

30. The display will read "DEL ALPHA = (SOME NUMBER) ." 

Note the n\amber and press "CONTINUE EXECUTE." 

31. The display will read "BETA 1 = (SOME NUMBER) ." Note 
the number and press "CONTINUE EXECUTE." 

32. The display will read "INPUT ZETA PO, FIG. 15, 1174 
VA^?" Turn to Fig. 15 in Ref. 16. Select the value of 
ZETA PO which corresponds to the values of DEL ALPHA, 
BETA 1; and ALPHA EXIT. press ("ZETA PO) EXECUTE." 

33. The display will read "RE ROTOR = (SOME NUMBER)." 

Note this number and press "CONTINUE EXECUTE." 

34. The display will read "INPUT K-RE, FIG. 23, GA1074VA2." 
Turn to Fig. 23 in Ref. 18 and select the value of K__ 

Khi 

that corresponds to the rotor blade Reynolds niimber 
given in step 33. Press "(K-RE) EXECUTE." 
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35. The display will read "M2 (IS) = (SOME NUMBER)." 

Note this mjmber and press "CONTINUE EXECUTE." 

36. The display will read "INPUT KM, FIG. 18, 1174VAj^." 

Turn to Fig. 18 in Ref. 16, and select the value of KM 
that corresponds to the blade exit Mach Number as 
calculated from isentropic conditions. Press "(KM) 
EXECUTE. " 

38. The program will CHAIN with "TTR 6." The display will 
read "ENTER NEXT RED REC. #." Input the same record 
number from step 12 and press "EXECUTE." 

39. The display will read "STORE RED. DATA ? YES = 1, 

NO = 0." Press "1 EXECUTE" if you desire reduced data 
to be stored and used for plotting and tabulation. 

Press "0 EXECUTE" if you do not desire to store reduced 
data. Go to step 41. 

40. You have elected to store the reduced data in the Mass 
Memory, and the words "REDUCED DATA IS STORED IN FILE/ 
REDDAT/RECORD (NUMBER FROM STEP 12)." will be printed. 

41. Electing not to store the reduced data the words "THIS 
DATA WAS NOT STORED" will be printed. 

42. The display will read " (CONT) (EXEC) MORE DATA - OR 
(END) ." 

43. If more data points are to be reduced, press "CONTINUE 
EXECUTE." The display will read "TAB RAW DATA ? YES = 1 
NO = 0." If you desire not to tabulate raw data, press 
"0 EXECUTE." Then the display will read "TAB REDUCED 
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DATA ? YES = 1, NO = 0.” If you desire not to tabulate 
reduced data, press "0 EXECUTE." By inputting zeros 
for both of these tabulation questions the program will 
display "REDUCE NEXT POINT," and will CHAIN with "TTR 
1." Go to step 12. Note; Tabulation is usually 
elected after all points are reduced for a run. 

44. To tabulate the raw data press "1 EXECUTE." The program 

will GET "TTR 7." The display will read "EXECUTE RECORD 
#'S: LOWEST, HIGHEST." Input your first and last 

record numbers and press "EXECUTE." The display will 
read "OMIT 4 RECORDS. ENTER #'S OR O'S." This option 
permits you to omit any points. Input the record 
numbers and press "EXECUTE." 

45. The display will read, read "REQUEST DATA CH # : 0, 1 

OR 10." Press "1 EXECUTE" will cause the 48 measurements 
to be tabulated under the heading; TTR INPUT DATA . 
Afterwards, the calculator will request the next channel 
nxamber. Press "10 EXECUTE" will cause all other channels 
to be tabulated and printed with the same heading. Next, 
press "0 EXECUTE." The display will read "TAB REDUCED 
DATA ? YES = 1, NO = 0." Press "1 EXECUTE" and the 
program will GET "TTR 8." Go to step 46. Press "0 
EXECUTE" and the program will GET "TTR 1." Go to step 
12 . 

46. By electing to tabulate the reduced data, the program 
will GET "TTR 8." The display will read "ENTER LOWEST, 
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HIGHEST REC. #.” Input the first and last record 
numbers, and press "EXECUTE.” 

47. The display will read "ENTER 1ST AND LAST PT. #." 

Input the two points, and press "EXECUTE.' 

48. The display will read "ENTER RUN #." Input the run 
number, and press "EXECUTE." 

49. The velocities, Mach Number, and losses will be tabu- 
lated and printed. The tabulation will stop to allow 
the page to be removed from the printer. Then press 
"CONTINUE EXECUTE." 

50. Tabulation and printing will resume until "DATA TABULATION 
COMPLETE" is printed. The display will read "CHAIN 

WITH TAPE READER. " There is a 3 second pause , and the 
program will GET "TTR 1." Go to step 12. 
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CHAN- PORT MEASUREMENT X SCALE = UNIT ASSIGNMENT 
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TABLE B-2 



DEFINITION OF VARIABLES AND EQUATIONS 
A Array 

AO Area of Control Volume (Fig. B-2, B-3) 

Al Area of Control Volume 

A3 Alpha 1 

A4 Alpha 2 

A5 Decision Variables 

A8 Rotor Shroud Pressure Integrated Over Area 

A9 Rotor Shroud Pressures Integrated Over Area 

BO EQ A (8) Ref. 5 

B1 B2 Eq. A (16) Ref. 5 

B2 Beta 1 

B3 Beta 2 

B4 Units (2*CP*g) 

B5 Rotor Shroud Pressures Integrated Over Area 

B6 Rotor Shroud Pressures Integrated Over Area 

B7 Rotor Shroud Pressures Integrated Over Area 

B8 Del Alpha 

B9 Mixing Loss Coefficient, Referred to Average Kinetic 

Energy, Eq. 71, Ref. 16 

C Array 

CO Gas constant (ft-lbf/lbm-R) 

Cl Specific Heat (BTU/lbm-R) 

C2 Nozzle Diameter (in.) 
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C3 

C4 

C5 

C6 



2 

Gravitational Constant (32.174 Ibm-f t/lbf-sec ) 
Pipe Diameter (upstream) 



Ratio (D /D ) 
noz p 

Conversion factor ; 



0.16384 

/2.036 



Cl 


Discharge 


Coefficient 


C8 


Thermal Expansion Coefficient 


C9 


Compressibility Coefficient 


D 


Array 




DO 


See Table 


B-1 


Dl 


See Table 


B-1 


D2 


See Table 


B-1 


D3 


See Table 


B-1 


D4 


See Table 


B-1 


D5 


See Table 


B-1 


D6 


See Table 


B-1 


D7 


See Table 


B-1 


D8 


See Table 


B-1 


D9 


See Table 


B-1 


EO 


Effective 


Total-to-Static Efficiency Eq. A (18) Ref 


El 


Next record Number 


E2 


Input Variables 


E3 


See TRble 


B-1 


E4 


Profile Loss Coefficient of Turbine At Design Point 


E5 


Profile Loss Coefficient (? ) Eq. 46, Ref. 16 

P 


E6 


Secondary 


Flow Loss (?) Eq. 63, Ref. 16 

o 


E7 


Mixing Loss Coefficient (^^) Eq. 74, Ref. 16 



5 
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I 



E8 

E9 

FO 

FI 

F2 

F3 

F4 

F5 

F6 

F7 

F9 

G 

GO 

G1 

G2 

G6 

G7 

G8 

G9 

H 

HO 

HI 

H2 

H5 



Stator Loss Theoretical (Total Blading Losses) 
Eq. 11 , Ref. 16 

Rotor Loss Theoretical (Total Blading Losses) 
Eq. 11 , Ref. 16 

Overall Axial Force 

Partial Sum of Overall Axial Force 



^cl^cl 

^a^tot 

Stator Axial Force 
Closure Plate Force 

V 

Ratio: (Thickness to Spacing) Eq. 73, Ref. 16 

Mixing Loss Coefficient, Referred to Average Kinetic 
Energy, Eq. 71, Ref. 16 

Correction Factor for Blade Thickness, Eq. 47, Ref. 16 

Array 

Y - 1 
Y 

Input Variable 

X 



Y - 1 

Decision Variable 
Decision Variable 

Secondary Flow Loss, Eq. 63, Ref. 16 

Thickness to Spacing, Eq. 73, Ref. 16 

Pressure at H\ib 

See Table B-1 

Horsepower 

Referred Horsepower 

Correction Factor For Effects of Reynolds Number 
'■'Ke' 
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H6 Correction Factor For Blade Thickness (K.) Eq. 47, 

Ref. 16 ^ 

H8 Density Appendix B, Ref. 5 

H9 Reynolds Nimber of Rotor, Table B-2, Ref. 5 

I First Point (Used in Tabulation Program) 

10 Ratio; Non-Dimensional Pressures At Various Tap 

Locations In Stator (Fig. 5) - DO in Table B-1 

II Ratio: Non-Dimensional Pressures at Various Tap 

Locations In Stator (Fig. 5) - D1 in Table B-1 

12 Ratio: Non-Dimensional Pressures At Various Tap 

Locations In Stator (Fig. 5) - D2 in Table B-1 

13 Ratio: Non-Dimensional Pressures At Various Tap 

Locations In Stator (Fig. 5) - D3 in Table B-T 

14 Ratio: Non-Dimensional Pressures At Various Tap 

Locations In Stator (Fig. 5) - D4 in Table B-1 

15 Ratio: Non-Dimensional Pressures At Various Tap 

Locations In Stator (Fig. 5) - D5 in Table B-1 

16 Ratio: Non-Dimensional Pressures At Various Tap 

Locations In Stator (Fig. 5) - D6 in Table B-1 

17 Ratio: Non-Dimensional Pressures At Various Tap 

Locations In Stator (Fig. 5) - D7 in Table B-1 

18 Ratio: Non-Dimensional Pressures At Various Tap 

Locations In Stator (Fig. 5) - D8 in Table B-1 

19 Ratio: Non-Dimensional Pressures At Various Tap 

Locations In Stator (Fig. 5) - D9 in Table B-1 

J Mach # At Station 9 In Stator 

JO Joules Constant (778.16 Ft-lbf/BTU) 

J2 See Table B-1 

J4 Input Variable 

J5 Profile Loss Coefficient Eq. 46, Ref. 16 

J6 Velocity at Station 9 in Stator 

J7 Non-Dimensional Stator Pressure 
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J8 Velocity at EX-TIP 

J9 Velocity at Station 5 in Stator 

K Fourth Degree Polynomial For K.E Factor (Fig. 4) 

KO Viscosity (Sutherland Formula) 

K1 Kinetic Energy Factor (G. E. Report) 

K2 Beta (G. E. Report) 

K3 K X B X X Pj^ (KxK2xk6xQ0) 

K4 Isen tropic Head Coefficient Eq. a (22) , Ref. 5 

K5 Slope of K. E. Curve (Newtonian Method) 

K6 Area of Labyrinth 

K7 K. E. Factor (Convergence) 

K8 Second Degree Polynomial K^(R^) For Nozzle (Fig. 5) 

K9 See Table B-1 

L Counter 



LO 


a^ + (Xu2 - U2) 


^ Eq. A 


(21) 


, Ref. 5 






LI 


Rotor Loss Coefficient Eq. 


A (21) , Ref. 


5 




L2 


Mach # At 
(Fig. 5) 


Various 


Pressure 


Tap 


Locations 


in 


stator 


L3 


Mach # At 
(Fig. 5) 


Various 


Pressure 


Tap 


Locations 


in 


Stator 


L4 


Mach # At 
(Fig. 5) 


Various 


Pressure 


Tap 


Locations 


in 


Stator 


L5 


Mach # At 
(Fig. 5) 


Various 


Pressure 


Tap 


Locations 


in 


Stator 


L6 


Mach # At 
(Fig. 5) 


Various 


Pressure 


Tap 


Locations 


in 


Stator 


L7 


Mach # At 
(Fig. 5) 


Various 


Pressure 


Tap 


Locations 


in 


Stator 


L8 


Mach # At 
(Fig. 5) 


Various 


Pressure 


Tap 


Locations 


in 


Stator 
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L9 Mach # At Various Pressure Tap Locations in Stator 

(Fig. 5) 

MO Referred Moment 

Ml Mach # At Station 1 

m 2 Axial Mach # At Station 1 

M3 Mach # At Station 2 

M4 Axial Mach # At Station 2 



M5 


See 


Table 


B-1 




M6 


See 


Table 


B-1 




M7 


Referred Stator Moment 


M8 


Mach # At 


Exit 


- Tip 


M9 


Mach # At 


SS - 


Aft 


NO 


See 


Table 


B-1 




N1 


See 


Table 


B-1 




N2 


See 


Table 


B-1 




N3 


See 


Table 


B-1 




N4 


See 


Table 


B-1 




N5 


See 


Table 


B-1 




N6 


See 


Table 


B-1 




N7 


See 


Table 


B-1 




N8 


See 


Table 


B-1 




N9 


See 


Table 


B-1 




OO 


Incidence 


Loss 


Coefficient 


01 


See 


Table 


B-1 





02 Velocities at Various Pressure Tap Locations In 
Stator (Fig. 5) 

03 Velocities at Various Pressure Tap Locations In 
Stator (Fig. 5) 
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04 


Velocities at Various Pressure 
Stator (Fig. 5) 


Tap 


Locations 


In 


05 


Velocities at Various Pressure 
Stator (Fig. 5) 


Tap 


Locations 


In 


06 


Velocities at Various Pressure 
Stator (Fig. 5) 


Tap 


Locations 


In 


07 


Velocities at Various Pressure 
Stator (Fig. 5) 


Tap 


Locations 


In 


08 


Velocities at Various Pressure 
Stator (Fig. 5) 


Tap 


Locations 


In 


09 


Velocities at Various Pressure 
Stator (Fig. 5) 


Tap 


Locations 


In 


PO 


Reference Pressure (29.92 in. 


Hg) 






Pi 


P-HOOD/P-LAB (RATIO) 








P2 


^ 0 (sin a^) 


(Alpha Input) 


P3 ■ 


pyp Eq. A (10) Ref. 5 




(K Input) 




P4 


w 

^2 








P5 


Pl/Pti 








P6 


Pressure Ratio (P. /P«) 

to ^ ^ 








P7 


See Table B-1 








P9 


Average Pressure At Rotor Exit 








QO 


See Table B-1 








Q1 


See Table B-1 








Q2 


^to^^atm 








Q3 


See Table B-1 








Q4 


See Table B-1 








Q5 


Decision Variable 








R 


Array 








Re 


Reynolds Number (Nozzle) 
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RO 

R1 

R2 

R3 

R4 

R8 

R9 

S 

S 

50 

51 

52 

53 

54 

55 

56 

57 

58 

59 
T 

TO 

T2 

T3 

T5 

T6 

T7 

T8 



Referred RPM (N//9) 

Reynolds Number (Labyrinth) 

R Eq. A (4) , Ref. 5 

Theoretical Degree of Reaction Eq. A (23) , 
Actual Degree of Reaction Eq. A (24) , Ref. 
See Table B-1 
Ratio (14.696/P^^) 

Array 

Temperature Subroutine Variable 
^1 

Temperature Value 
Temperature Value 
See Table B-1 
See Tcible B-1 
See Table B-1 
See Table B-1 
See Table B-1 
See Table B-1 
See Table B-1 
Pressure at Tip 

Reference Temperature (518.7 ®R) 

T^^ (Upstream Temperature ®R) 

/F Temperature Ratio 

Ttn (Nozzle Temperature °R) 

T^p (Plenum Temperature ®R) 

Tg (Fig. B-4) 

T^ (Fig. B-4) 



Ref. 5 
5 
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T9 


(Pig. B-4) 




U 


Array 




UO 


Mixing Loss Coefficient, Eq. 74 


, Ref. 16 


U1 


(Dq + D^)/48 




U2 


(Dq + D^)/48 




U3 


Dimensionless Velocity 




U4 


Dimensionless Velocity 




U5 


Factor Cf P* 37, Ref. 16 




U6 


Loss Due to Tip Clearance Flow, 


Eq. 75, Ref. 16 


U7 


Tip Clearance Loss Coefficient, 


Referred To Isen 




tropic Conditions Eq. 76, Ref. 


16 


U8 






U9 


stator Flow Function 


■ 


V 


Array 




VO 


Tangential Rotor Velocity 




VI 


Limiting Velocity (V^^) 




V2 


V ^2 Average Tangential Velocity 


Eq. A (11) , Ref. 


V3 


1 — 1 
nJ 

> 




V4 


CM 




V5 


^2 




V6 


^1 




V7 


Eq. A (10) , Ref. 5 




V8 


Stator Blockage Factor 




W 


Array 




WO 


Total Flow Rate (W = W “ W, , ) 

noz lab 


W1 


Labyrinth Flow Rate Based on K. 
Report) 


E. Factor (G. E. 
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W2 Labyrinth Convergence Flow Rate 

W3 Flow Rate (ASME Equation) 

W4 Ratio (WO/g) 

W5 Referred Flow Rate 

W6 (See Fig. B-5) 

W7 W 2 (See Fig. B-5) 

W8 W*^ = Cos (Bj^ - 62) 

W9 W ^2 ^2ee Fig. B-5) 

X Array 

XO Dimensionless Velocity X^^^ 

XI - Eq. A (7) , Ref. 5 (K Input) 

XI X^j^ (Alpha Input) 

X2 X=X^ + X^ 

u2 a 

X3 1/A PdA 

X4 X^2 ^<5* ^ 5 

X5 X^ Eq. A (15) , Ref. 5 

X6 X 2 

X7 ^/^to^ (Alpha Input) 

X8 y(2K-1) + 1 

X9 Denominator of Eq. 63, Ref. 16 

Y Array 



YO 

Y1 

Y2 




P ^ 
-Cp^) Y 

to 




izi 
) Y 




For 10 



For II 



For 12 
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For 13 



For 14 



For 15 



For 16 



For 17 



For 18 



For 19 



z 


Array 




zo 


Cl Eq. A (8a) , Ref. 5 




Z1 


Stator Loss Coefficient 




Z2 


C2 Eq. A (16) , Ref. 5 




Z3 


Rotor Loss Coefficient Eq. 


A (21) , Ref. 


Z6 


Non-Dimensional Pressure At 


Exit - Hub 


Z7 


Loss Due Tip Clearance, Eq. 


75, Ref. 16 


Z8 


Rotor Loss Coefficient (Other Factors) 
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PRDGRHM • TABLE B-3 




SCORE SHEET: RECORD OF VRRIHBLES USED. 
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FIG. B-1 DATA REDUCTION SCHEMATIC 
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FIG. B-1 DATA REDUCTION SCHEMATIC 



TTR 5 

1.) PRINT REDUCED DATA 



ik 

TTR 6 

1.) STORE REDUCED DATA 



TTR 7 

1.) TABULATION OF RAW DATA 



iU I 

TTR 8 

1.) TABULATION OF REDUCED DATi. 
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FIGURE B-3 CONTROL VOLUME b 
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FIGURE B-4 THERMODYNAMIC PROCESS OF FLUID 
IN AN AXIAL TURBINE STAGE 
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FIGURE 3-5 VELOCITY DIAGRAM 
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